Statistical analysis of coherent structures in transitional pipe flow 
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Numerical and experimental studies of transitional pipe flow have shown the prevalence of coherent 
flow structures that are dominated by downstream vortices. They attract special attention because 
they contribute predominantly to the increase of the Reynolds stresses in turbulent flow. In the 
present study we introduce a convenient detector for these coherent states, calculate the fraction of 
time the structures appear in the flow, and present a Markov model for the transition between the 
structures. The fraction of states that show vortical structures exceeds 24% for a Reynolds number 
of about TZe — 2200, and it decreases to about 20% for TZe — 2500. The Markov model for the 
transition between these states is in good agreement with the observed fraction of states, and in 
reasonable agreement with the prediction for their persistence. It provides insight into dominant 
qualitative changes of the flow when increasing the Reynolds number. 

PACS numbers: 47.27.Cn, 47.27.eb, 47.27.nf, 47.27.De 
Keywords: pipe flow, coherent structures, Markov model 



I. INTRODUCTION 

The visualization of turbulent flows and boundary lay- 
ers via sophisticated experimental methods Hke particle 
imaging velocimetry has led to the identification of a 
rich variety of prominent coherent structures, such as 
waves .streaks, hairpin vortices and lambda vortices[ll, 
0, [^, These extended coherent structures signifi- 

cantly influence large scale momentum transport and 
hence Reynolds stresses. As a consequence, they figure 
prominently in turbulence research. 

Studies on internal flows in confined geometries have 
highlighted the dominant role of structures containing 
pronounced downstream vortices and have led to the pro- 
posal of a three-step self re gene rating mechanism for tur- 
bulence @, i, 0, H, M, d [ntlil, m • Downstream vortices 
transport liquid across the mean shear gradient and cre- 
ate regions of fast or slow moving fluid, so-called high- 
and low-speed streaks. The streaks generated by this lift- 
up process become unstable to the formation of normal 
vortices, and through a nonlinear interaction mechanism 
the latter feed their energy back to downstream vortices. 
This closed self-regeneration mechanism appears to be a 
generic dynamical feature of turbulent shear flows. The 
process was identified in direct numerical simulations of 
plane Couette flow in narrow cells where transverse mod- 
ulations are constrained [1, d, [l^, but it can also be 
detected in time-correlation functions in fully turbulent 
flows 'l4\. 

In its purest form this self-regenerating cycle gives rise 
to a periodic solution to the equations of motion. How- 
ever, in most coherent structures the flow is not strictly 
periodic and always perturbed by background fluctu- 
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ations. Examples of exact coherent states have been 
given in simple models, where they correspond to peri- 
odic orbits 15, 16, l3) Ell, and, in the full flow, through 
the numerical identification of exact coherent states in 
channel fiows [H IH, M, HE and 

travelling waves in pipe flow [la [li [IJ In all these 
cases, the coherent structures are dominated by pairs 
of counter-rotating downstream vortices and associated 
streaks which are regularly arranged in azimuthal direc- 
tion. The flow fields are invariant under discrete rota- 
tions around the pipe axis. 

Since all exact coherent states constructed so far are 
linearly unstable it came as a surprise that they could be 
directly observed in experiments [30| ■ In this work we fol- 
low up on this experimental observation with a study of 
the appearance and persistence of these structures in nu- 
merical simulations of pipe fiow. In particular, we show 
how they can be detected, how frequently they appear, 
and how long they persist. 

The travelling waves observed in pipe flow are of par- 
ticular interest because they are believed to form a back- 
bone for the turbulent dynamics near the onset of tur- 
bulence. Since the laminar profile in these flows is lin- 
early stable for all Reynolds numbers 0, [H, [H, [H, [H, 
m, llSj nil the transition cannot proceed through states 
bifurcating from the laminar profile. The turbulent mo- 
tion which in many pipe-flow experiments is observed for 
Reynolds numbers beyond about 2000, must hence arise 
via a nonlinear transition scenario 0, [§, [3§, Hi], \^ . 
The travelling waves are then the simplest persistent non- 
linear structures around which the turbulent dynamics 
can form. Together with their stable and unstable man- 
ifolds they can give rise to the basic building blocks of 
chaotic dynamics, such as hyperbolic tangles and Smale 
horseshoes. While it is unlikely that one will be able to 
identify an individual travelling wave in a time series, it 
is possible to identify a visit to their neighborhood, as 
identified by the appearance of similar structures in the 
flow. 
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In the present paper we propose a way to detect the 
visits to the neighborhoods of coherent states, and use it 
to infer information about the structures underlying tur- 
bulence. To distinguish different parts of state space and 
different flow topologies, we introduce projections onto 
lower dimensional subspaces that capture salient features 
of classes of coherent states, and study the recurrences 
to these subspaces: This is weaker than identifying in- 
dividual travelling waves but sufficient to discriminate 
between various flow regimes. On the technical side, the 
reduction in resolution also lowers the requirements on 
the length of the time traces and helps to improve the 
statistical significance. 

The outline of the paper is as follows. In Sect. [Til we 
briefly describe the spectral code underlying the simu- 
lation of the flow. In Sect. IIIII we describe the projec- 
tion used to detect and characterize the coherent struc- 
tures in direct numerical simulations of pipe flow close to 
the threshold of turbulence. In Sect. IIVI we analyst the 
statistics of the occurrence of coherent structures, and in 
Sect. |V] we explore their physical properties. The paper 
closes with a discussion and outlook in Sect. I VII 



II. SIMULATION OF PIPE FLOW 

We consider an incompressible Newtonian liquid in a 
pipe of circular cross section subject to no-slip boundary 
conditions at the walls. The flow is forced by a uniform 
pressure gradient which is adjusted to keep the flux con- 
stant at any instant of time [30l [43l [4^ . In other words 
the integrated volume flux through a cross section of the 
pipe is constant, and the Reynolds number 



Tie 



(1) 



is externally controlled in order to be independent of the 
flow state of the liquid. Here {u^) denotes the mean 
downstream velocity, R is the pipe radius and v the kine- 
matic viscosity of the liquid. In our simulations the pipe 
is L = lOi? long, and we use periodic boundary condi- 
tions in the downstream direction: physically, this corre- 
sponds to a numerical representation of the interior of a 
turbulent patch. 

The Navier-Stokes equations are written in cylindrical 
coordinates (r, if, z) and solved with a pseudo-spectral 
scheme. In doing so dimcnsionless units where lengths 
are measured in units of the radius of the pipe and ve- 
locities in units of twice the mean downstream velocity 
{i.e., the centerline velocity of the equivalent parabolic 
laminar profile) are used. Time is measured in units of 
i?/2K). 

All three components of the velocity field (ur,u^,Uz) 
are decomposed into Fourier modes in azimuthal and 
downstream direction. Chebyshev polynomials are used 
for expansion in radial direction. The velocity field is 




FIG. 1: The coherent vortices are expected to be most promi- 
nent in cross sections perpendicular to the pipe axis, as the 
shaded plane in this figure. In experiments, the velocity fields 
in this cross section are obtained by steroscopic particle image 
velocimetry [SCt] . 



thus written as 




Here the spectral basis functions are 
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wher e T, - denotes the j*'' normalized Chebyshev polyno- 
mial [U nil , and kz = In physical collocation space 
the velocity fields are represented by the values of the 
fields at the corresponding Gauss-Lobatto grid points. 

A 4th-5th-order Runge-Kutta-Fehlberg scheme with 
adaptive step-size control is used to evolve the solution 
in time [i^l, and the action of the Navier-Stokes operator 
is computed via a pseudo-spectral scheme. The transfor- 
mation between spectral and physical space required by 
the pseudo-spectral scheme is performed by FFT-based 
routines. Constraints (incompressibility, regularity and 
analyticity) as well as no-slip boundary conditions are 
enforced by a Lagrangian projection mechanism [48[ . 

The simulations presented in this work are carried out 
with n Fourier modes in azimuthal and m Fourier modes 
in downstream direction, where -I- < 1 . Conse- 
quently, we consider up to 49 Fourier modes in azimuthal 
and up to 45 in the downstream direction. 47 Cheby- 
shev polynomials are used for the expansion in the radial 
direction, adding up to 3 x 49 x 23 x 47 Ri 1.6 x 10^ 
components. 



III. DETECTION OF COHERENT 
STRUCTURES 

The travelling waves [l^, we want to detect are 
dominated by vortices aligned along the axis, and corre- 
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FIG. 2: (color online) Deviation u~ (u)t of the instantaneous 
velocity field u from the mean turbulent profile (u)t for a 
pipe flow at TZe — 2200. The colors indicate the downstream 
velocity component according to the scale specifled by the 
color bar to the right, and the in-plane velocity components 
are indicated by arrows. The two panels show (a) a case where 
no clear structure is observed, and (b) one with a four-fold 
streak. 



spending streaks in the downstream velocity components. 
The downstream vortices and streaks are most prominent 
in cross sections of the pipe perpendicular to the axis. As 
in the experiments [33|, where stereoscopic particle im- 
age velocimetry was used to extract the velocity fields, 
we will focus on the velocity fields in cross sections per- 
pendicular to the pipe axis (Fig. [T|). For the travelling 
waves it makes no difference whether we focus on one 
cross section and follow the time evolution or whether we 
freeze the flow at one instance of time and move the cross 
section along the axis. The same applies for a transient 
appearance of these structures: in a fixed cross section 
they will come and go, and in a frozen flow they would 
be present in some regions along the axis and absent in 
others. In the analysis presented below we work, as in 
the experiments, with the time evolution in cross sections 
at a fixed position in the lab frame. Typical examples of 
cross sections with high- and low-speed streaks, i.e., of 
regions of high and low downstream velocity, are shown 
in Fig. [21 The structures are best visible when a reference 
profile is subtracted. In previous works (s^l the laminar 
profile with equal mean velocity was subtracted. Here we 
use the mean turbulent profile. It is obtained as the av- 
erage over azimuthal angle and time of the downstream 
velocity at a fixed radius. 



A. Characterizing the symmetry of coherent states 

As mentioned in the introduction the coherent trav- 
elling waves identified so far have highly symmetric ar- 
rangements of vortex pairs. By transporting fast liquid 
from the center to the walls and slow liquid from the 
wall to the center region, these pairs of vortices gener- 
ate elongated regions of fast and slow moving liquid. We 
therefore focus on the appearance of symmetric arrange- 
ments of high- and low-speed streaks schematically in- 



FIG. 3: (color online) Sketch of the regular arrangement of 
high- (dark red) and low-speed (light blue) streaks in coherent 
structures. When analyzed at a fixed radial position close 
to the wall (green dashed line at radius 0.81), all currently 
known travelling-wave solutions show high-speed streaks that 
are equidistantly arranged on the circumference, i.e., they 
show an A'^-fold rotational symmetry. Typical states contain 
N low speed streaks close to the center, and 2N high speed 
streaks close to the wall f27|. However, also states containing 
A'' low and A*' high speed streaks were found [29l |. 



dicated in Fig. [31 The travelling waves also show that 
the high-speed streaks close to the walls are fairly stable 
and do not move much in the azimuthal direction over 
one period. This simplifies their detection amidst the 
fluctuations of the total velocity field. 

The rotational symmetry of the pipe entails that pat- 
terns should be considered identical when they only differ 
by a global rotation around the pipe axis. A detector for 
coherent states should take this into account and be in- 
variant under global rotations. We therefore suggest to 
use an azimuthal correlation of the downstream velocity 
Uz at a chosen radius r and axial position z. 



C{(l)) = {uz (r, ip, z) Uz (r, + 0, z))^p 



(4) 



where (•)(^ denotes averaging over ip. By a straightfor- 
ward calculation one verifies that this correlation func- 
tion is invariant under global rotations. Moreover, it re- 
liably uncovers periodic structures in the azimuthal di- 
rection. 

Whenever the system approaches a coherent state 
showing N high-speed streaks close to the wall, the cor- 
relation function C((/)) shows N peaks separated by an 
angular displacement 2n/N . In particular, the four-fold 
structure of the downstream velocity field. Fig. [2b results 
in a clear four-fold structure of the correlation function, 
which is shown in Fig. [3^. In addition to the autocor- 
relation peak at = the correlation function shows 
peaks at = ±-|,7r. Similarly, a flow with a three-fold 
symmetry gives rise to peaks at (/> = 0, ±^ (cf. Fig. [Ha). 

By following C{(j>) in time one can detect the lifetimes 
of structures, their decay, and the subsequent emergence 
of new patterns. An example is given in Fig. [5l which 
shows the decay of a four-streak state and the emergence 
of a six-streak state within about 1 pipe radius. 
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FIG. 4: (color online) Azimuthal correlation functions evalu- 
ated at r = 0.81 for the velocity fields shown in Fig. [5] When 
no clear structure is observed in the cross section, the correla- 
tion function only shows an autocorrelation peak at = (red 
circles), (a) When a four-fold symmetry is present, e.g. for 
t = 1217.2, the correlation functions has additional peaks at 
= ±-| and 77 (yellow boxes), (b) For a three- fold symmetric 
state, the additional peaks appear at = 



B. Automated structure detection 

The correlation function C{(f>) signals the proximity of 
the flow to a coherent state by evenly spaced peaks. Its 
derivatives highlight both minima and maxima of the cor- 
relation function (see Fig. [5]) and emphasize flow struc- 
tures of comparable (azimuthal) streak gradients. Since 
C(0) is an even function in cj), its derivative is odd. It 
should have an substantial overlap with the sine-function 
of the appropriate periodicity. In order to automatically 
detect evenly spaced maxima and in order to count their 
number we therefore define the scalar measures via 
a scalar product of the derivative of the correlation func- 
tion and sin(A^(/)) [sl ]. 




ZNit) 



d^C{(j)) sin(iV<^) dcj) 



(5) 



0,01 



o.oo- 



■0,01 



FIG. 5; (color online) Azimuthal correlation function plotted 
as a function of downstream position in the pipe. One clearly 
observes the transition form a four-streak state to a six-streak 
state. The transition is quite sharp and happens within a 
spatial range of a single pipe radius. 



This reduction of information to a scalar quantity con- 
tains one parameter, the radius r at which the correla- 
tion functions are determined. For the Reynolds numbers 
considered here we find r = 0.81 to be convenient. At 
this radius, which is indicated by a dashed green line in 
Fig. [3l the coherent structures under investigation show 
a pronounced regular arrangement of high-speed streaks. 

By following time traces oi Z^q for different TV we can 
study the prevalence of structures of certain multiplicity 
and the transitions between them. Examples are given 
in Fig. [6l The top frames show d^C{(l)), the derivative 
of the azimuthal correlator with respect to the angular 
coordinate (f), as a function of the azimuthal coordinate 4> 
and the time t. The four-fold structures have eight zeros 
in their derivative (from four maxima and four minima), 
and the six- fold structures have twelve zeros. Parallel 
nodal lines indicate the presence of these structures for 
times of about ten natural time units. 

The lower frames in Fig. [6] show the time evolution 
of the corresponding scalar projectors Zjv- The indica- 
tor Z4 shows pronounced peaks when the four-fold sym- 
metric patterns are observed in the correlation function 
and Zq peaks when the six-fold structures appear; con- 
versely, one is small when the other one is large. One 
also notes considerable fiuctuations due to the residual 
background turbulence. In general, values of Zn smaller 
than about 0.01 cannot be considered significant indica- 
tors of a structure and belong to background fluctuations. 
On the other hand, comparison of the top and bottom 
frames in Fig. [6] suggests that a threshold Zn > 0.013 
signifles the presence of coherent structures with iV-fold 
symmetry. 

Armed with this threshold, we collapse the scalar time 
series Z^it) for iV = 2, . . . , 8 to a single discrete indi- 
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FIG. 6: (color online) The derivative of the correlation function dii,C{(j)) as a function of (f> and time t (top), and of the 
corresponding scalar measures Z4 (blue solid line) and Ze (red dashed line) in bottom part. The color coding in the top graphs 
runs linearly from —0.005 (blue) to 0.005 (red). Nodal lines appear in white. Besides irregular, featureless correlation functions 
at e.g. t = 1200 . . . 1208 and around t — 1230, there are long stretches of time where the function shows a distinct four-fold 
(e.g. a.t t — 1210 . . . 1220) and six-fold symmetry (e.g. around t = 1334 and 1350), respectively. 



cator, N{t), which assigns to a cross section at time t 
the number N of symmetric streaks and corresponding 
vortices it contains. N takes the values 0, 2, 3, . . . , 8, 
where = is assigned to cases where all Zjv remain 
below the threshold. The maximal value 8 is an empiri- 
cal limit, in that states with eight or more vortices were 
rarely realized. 



IV. STATISTICAL ANALYSIS OF THE TIME 
SERIES 

Based on the time series Z^it) we now explore the 
statistical properties of the occurrence of coherent struc- 
tures in pipe flow. The aim of this statistical analysis is 
twofold: we want to see how frequently structures of a 
certain multiplicity are present and we want to study the 
extend to which a Markov approximation can describe 
the switching between states. 



A. Probability distribution of coherent states 

Fig. [7] shows the probabilities of detecting a coherent 
state of A^-fold symmetry in time series taken at different 
Reynolds numbers TZe close to the transition to turbu- 
lence. For TZe — 2200 about 24% of all cross sections fall 
into the categories = 3, 4, 5, and 6. For TZe — 2500, 
the fraction decreases slightly to about 20%. This high 
fraction explains the ease with which coherent structures 
were picked out of experimental cross sections f3(|, and 
underlines their significance as building blocks of the tur- 
bulence in the transition region. 

With increasing Reynolds number the weight of states 
with large N increases. These structures are much closer 
to the walls where they give rise to steeper gradients in 
radial and azimuthal direction and consequently larger 
friction. As these structures have more spatial degrees of 
freedom, it is less likely that they appear in perfect sym- 
metry. Hence, their correlators have smaller amplitudes, 
and it would be interesting in forthcoming work to probe 
for the structures with a localized correlator. 
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FIG. 7: (color online) Comparison of the statistical weight of coherent structures with A'^-fold symmetry in transiently chaotic 
time series of flows with different Reynolds numbers: (a) Tie — 2200, (b) TZe = 2350, (c) TZe = 2500. The red bars in the front 
are directly calculated from the simulation output, and the green bars in the background are the prediction of the Markov 
model. 



(a) (b) (c) 




FIG. 8: (color online) Distribution of the persistence time of coherent structures with A''-fold symmetry for the data also shown 
in Fig. [7] Following Eq. (|9]) the slope of the straight lines is determined from the diagonal elements of the transfer matrix. 
Time is measured in units of the sampling time r = 1.4. (A'^ = 4: open blue circles and blue dotted line, N = 5: red filled 
boxes and solid red line, A'' = 6: open green triangles and green dashed line.) 



B. Markov model for transitions 

The typical persistence time of a pattern in Fig. [S] is 
about 5 to 10 time units, and the transition between the 
four-streak and six-streak state shown in Fig. [5] takes 
about one time unit. When discretizing time in order to 
describe the transitions between different patterns, the 
sampling time scale should therefore not be much longer 
than about 5. Otherwise one misses states. On the other 
hand, if the time steps are much shorter than unity, one 
begins to probe the continuity of the time evolution. As 
representative examples in this interval we explored the 
discrete dynamics of discretized time sequences with a 



time spacing of t = 1.4 and of r = 2.4. Since different t 
lead to results which cannot be distinguished within our 
error margins, we will in the following present data for 
T — lA only. 

By considering the underlying flow at multiples of the 
time unit t its continuous dynamics is transformed into 
a discrete time-series. The conditional probability that 
one encounters an A^'-streak state in the following snap- 
shot, when currently facing an TV-streak state defines a 
transition matrix T^i n. Its indices N and A'^' take the 
values (when there is no streak), and A^ = 2 . . . 8 when 
Zjv exceeds its threshold value. For the Reynolds number 
TZe = 2200 we find 
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/0.90 0.26 

0.00 0.73 

0.02 0.00 

0.04 0.01 

0.03 0.00 

0.01 0.00 

0.00 0.00 

Vo.OO 0.00 



0.26 0.27 

0.00 0.00 

0.72 0.01 

0.01 0.71 

0.01 0.01 

0.00 0.00 

0.00 0.00 

0.00 0.00 



0.38 0.55 

0.00 0.00 

0.01 0.01 

0.03 0.03 

0.57 0.03 

0.01 0.38 

0.00 0.00 

0.00 0.00 



0.71 1.00 \ 

0.00 0.00 

0.02 0.00 

0.07 0.00 

0.00 0.00 

0.04 0.00 

0.16 0.00 

0.00 0.00 / 



The columns of the matrices add up to 1 because each 
state has to go to one of the eight admissible states in 
the next time step, 

5] Tn'n = 1. (7) 

Af'=0,2,3...8 

More than 10 000 independent snapshots (17 202 for 
TZe = 2200, and more than 15 000 snapshots at the higher 
TZe = 2350 and TZe — 2500) were analyzed in our statis- 
tics. For the lowest Reynolds number the = 8 class is 
observed only a single time, and it immediately relaxed 
into the iV = state (cf. rightmost column of Tij). De- 
spite its rare occurrence, the N = 8 state is included in 
the analysis because its statistical weight increases with 
Reynolds number: It reaches 0.4% at TZe — 2500. 

C. Invariant distribution and life time of coherent 
states 

To check that the Markovian dynamics generated by 
the transition matrices faithfully represents the continu- 
ous dynamics, we first calculate the invariant probability 
distribution e, defined as the eigenvector to the eigen- 
value 1, i.e. 

e = Te . (8) 

Fig. [7| shows that the cn faithfully reproduce the relative 
frequencies in the original data. A comparison of the 
histograms for the different TZe shows that the number 
of visited coherent states and the complexity of the re- 
sulting flow patterns increase when the Reynolds number 
increases. 

Except for = 8 the highest transfer probabilities 
in each column appear along the diagonal of T. These 
elements describe the persistence of flow patters from 
time step to time step. Therefore, the probability density 
function p{N, n) to observe the pattern for n consecutive 
time steps scales like 

p{N,n) ^ {TNNr . (9) 

Fig. [5] shows data for the life time calculated from di- 
rect numerical simulation of the flow, together with 
the prediction from the Markov model, which is shown 



as straight lines in the semilogarithmic plot of life- 
times. Since long persistence times are exponentially 
suppressed this comparison requires very long time se- 
ries to check the prediction with reasonable statistical 
accuracy. Within these limitations there is a very good 
agreement between the data and the prediction. 

V. PHYSICAL PROPERTIES OF DETECTED 
STATES 

The different flow patterns also affect the different ve- 
locity and fluctuation statistics. As examples we con- 
sider the Reynolds stresses s^z = {uzU^}, Srr = (urUr) 
and Szr = (uzUr). Taking averages over 0, but not over 
time, provides probability distribution function (pdf) of 
temporal variations of these quantities (dashed lines in 
Fig. [5] tagged as combined), as well as conditional pdfs 
referring to states with a fixed number of streaks (solid 
lines tagged as state 3. . . state 6) and the turbulent un- 
structured state (solid lines tagged as state 0) . The over- 
all pdf can thus be decomposed into contributions of the 
previously discussed high-symmetry coherent states and 
a turbulent remainder (state 0). To emphasize the role 
played by the coherent states in changing the shape of the 
distribution of the considered component of the Reynolds 
stress the abscissa is always normalized to its overall tem- 
poral average. For instance, Szz is normalized by its av- 
erage si7, and the resulting normalized stress is denoted 
Szz = Szz/ Szz- By definition the mean of the Szz distribu- 
tion is therefore unity. However, the conditional pdfs for 
specific states will in general have means different from 
one. If the mean is larger than one, the state shows — 
on average — larger stress components than the tempo- 
ral average value of the component. Table [J lists both 
the absolute and the normalized mean values of all pdfs 
shown in Fig. [SI 

A. PDFs at fixed TZe 

From a physical point of view the interest of the de- 
compositions of the total pdf into conditional ones for 
turbulent and individual coherent states lies in the in- 
sight it gives into how the coherent states contribute to 
the exceptional statistics of fiuctuations in turbulent flow. 
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FIG. 9: (color online) The axial velocity fluctuations (top) , radial velocity fluctuations (middle) , and radial momentum transport 
(bottom) for flows at Reynolds numbers TZe — 2200 (left), Tie — 2350 (center) and TZe — 2500 (right). Different lines refer 
to the the overall time-averaged signal (dashed black, left axis), and the one averaged only over turbulent states, where no 
streaks are detected (solid black, left axes). The lines in the other colors give the respective contributions of states with a 
given number of streaks (right axis). The histograms are normalized with respect to their integral, i.e., the overall distribution 
(dashed line) is normalized to unity, and all other distributions (given by solid lines) add up to the overall distribution. Their 
norm consequently amounts to the weights shown in the histograms in Fig. [T] and the position of their maxima indicate for 
which values they most strongly contribute to the overall signal. 



We first consider the decomposition of the pdfs at fixed 
TZe, i.e., we discuss the trends in the mean of the data 
shown in individual panels of Fig. [5) 

On the average the detected coherent states generate 
much stronger Reynolds stresses than those found for the 
unstructured turbulent state 0. Consequently the coher- 
ent states shift the means and the maxima of the com- 
bined pdf to slightly larger values. Compared to the pdf 
of state (dashed black line) the coherent states add a 
fat tail to the combined pdf (black solid line) on the side 



of larger values of the stresses. In order to gain insight 
into the mutual importance of the different states we dis- 
cuss the trends in the mean and maxima as a function of 
the number of streaks N. 

The normalized stress Szz characterizes the intensity 
of streak structures in the flow field by estimating their 
downstream velocity. The maxima and mean values of 
its pdf decrease in the order A'' = 3, 4, 5 and 6. This can 
be interpreted as follows: The product of typical gradi- 
ents of Uz with the length scale over which the gradients 
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TABLE I; The temporal mean (tot) of the Reynolds stresses 

Szz ~ {UzUz), Srr ~ {UrUr) , and Szr ~ {UzUr) in UnitS of 

4(uz)^, and those of the corresponding conditional pdfs for 
disordered motion (state 0) and coherent states with A'^ — 
3 ... 6 streaks, respectively. In addition also the corresponding 
relative values Sij are given which are normalized with respect 
to the overall temporal mean of the considered component of 
the Reynolds stress. The related pdfs are shown in Fig. [9] 




FIG. 10: (color online) The solid black line shows the contri- 
bution of the states with steaks to the shift of the expectation 
value of Szr when TZe is increased from 2200 to 2350 and from 
2350 to 2500, respectively. The shift results from those in the 
restricted pdfs for states with 3, 4, 5 and 6 streaks, which 
are shown by broken lines with colors matching the choice 
adopted in Fig. [5] 



decrease in the position of the maxima in the sequence 
iV = 6, 5, 4 and 3. 



B. Drift of the mean with TZe 



persist is of the order of magnitude of the typical velocity 
fluctuation in downstream direction. Consequently, the 
azimuthal components of the gradients of Uz are of the 
same order of magnitude in all coherent states, and their 
typical length scale decreases like N~^. 

The radial component §„ measures the typical fluctu- 
ations of the radial velocity component, i.e., it charac- 
terizes the strength of the vortices. For this stress there 
also is a clear trend in the position of the maxima and 
mean values with N, but with the sequence reversed: 
the highest value for the maximum appears for TV = 6, 
and it decreases towards A'' = 3. This finding suggests 
that stronger vortices are needed to maintain the smaller 
streaks in coherent states with larger A^. 

From a physical point of view the Reynolds stress Szr is 
the most interesting of the three quantities. After all, it 
reflects the strength of the radial momentum transport. 
Hence it provides direct insight in the friction factor in 
the turbulent flow '49*1, and it also immediately reflects 
the role of the coherent states in the flattening of the 
laminar flow proflle in radial direction. In view of the op- 
posite scaling of the radial and axial velocity components 
observed in Srr and Szz, respectively, its A^ dependence 
results from a most subtle balance. Indeed, the coun- 
teracting trends almost cancel, leaving only a very weak 



In order to explore how the components of the 
Reynolds stress change with Reynolds number, and 
which physical effects generate the observed trends, one 
observes that the mean a; of a combined pdf P{x) = 
J2n ^nPn{x) with J dx P{x) = 1, / dxPN^x) = 1, and 
X^AT Cat = 1 is the weighted average of the means xn of 
the conditional distributions Pn{x), 

X— dx xP{x) =''^^eN / dxxPN{x) = cnXn 

J N •' N 

In Fig. [HI the conditional pdfs eMPwix) are plotted to- 
gether with their sum P{x) for x = Szz, Srr, and Szr, 
respectively, and the abscissa is scaled such that x = 1. 
The shift in the mean of P{x) therefore arises as an av- 
erage of the distance of the mean xn from unity with 
weights cn previously discussed in the framework of the 
Markov model (cf. Fig. [7]) . There are two physical effects 
underlying the observed changes in the statistics with 
TZe: (1) the change of the mean of conditional pdfs of 
the different states, and (2) the change in the statistical 
weights of the states. We will disentangle these contri- 
butions now for the physically most interesting case of 

Both visual inspection of the conditional pdf in Fig. [H] 
(bottom row), and the values of the normalized mean val- 
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ues in Table |T] show that there only is a slight drift of the 
coherent states' pdfs with TZe. In contrast, as observed 
upon discussing Fig. [7] their weights show pronounced 
changes. The nontrivial evolution of their weights with 
TZe suggests that the coherent states contribute to the 
change of the overall mean mainly by the change of their 
statistical weights e^r. This becomes particularly clear 
when plotting the relative change Aszr/szr of the posi- 
tion of the mean when increasing TZe from 2200 to 2350 
and from 2350 to 2500, respectively (Fig.[Tni). In the first 
interval this change is dominated by the one of state 3 
while state 6 hardly contributes, and in the latter interval 
these two states take just the opposite roles. 

We thus conclude that our statistical analysis allows us 
to identify the contributions of individual coherent states 
to the anomalous statistics of turbulent pipe flow, and 
to disentangle the changes with TZe into changes of the 
statistical weights of the states, and the comparatively 
smaller ones due to the T^e-dependence of the properties 
of individual states. One can interpret these findings 
as a hint that turbulent transients close to TZe > 2000 
are dominantly influenced by coherent states with only 
few streaks. In contrast, at higher TZe successively more 
coherent states with larger number of streaks affect the 
time series. 



VI. DISCUSSION 

In this section we want to summarize the results from 
the present simulation of pipe flow and point to the parts 
that could be useful in analyzing other shear flows as well. 

The automatic detection algorithm for coherent states, 
which was first used in [30j and was expanded on here, is 
fairly robust. It can be generalized to other flows as well. 
The algorithm systematically searches for structures that 
show a symmetric azimuthal arrangement of high-speed 
streaks along the wall which is topologically very simi- 
lar to the one observed in exact coherent states reported 
in [23, HO]. The detection is based on a Fourier- 
mode decomposition of the radial velocity. Because the 
detected states have the same symmetry structure also 
in the other components of the velocity field, the results 
should be robust against details of the implementation of 
the detector. Different projectors constructed along the 
line outlined in Sect. IIIll should lead to identical results. 
For extensions to larger TZe it might, however, be valu- 
able to consider extensions to asymmetric exp ansions of 
the flow fleld, e.g. by using wavelet [13, [Elf or Gabor 
[53 , [53} representations to extract basic units of coherent 
structures which contain only a single pair of vortices. 

In principle, one can obtain more highly resolved and 
more accurate information about the statistical proper- 
ties of the flow by including more degrees of freedom 
and subsequent extensions of the subspace of projection. 
In practice, however, the refinements are limited by the 
available data set, because more degrees of freedom re- 
quire many more data points in order to guarantee sta- 



tistically reliable results. 

Irrespective of the chosen detector, the present method 
can be used to quantitatively analyze coherent structures 
in turbulent flow: the automatic projectors give infor- 
mation about the probability to observe certain coherent 
structures, their lifetimes, and the transitions between 
these different states. A number of observations can thus 
be made: 

1. Despite being unstable the coherent states carry 
a considerable statistical weight of ~ 20%. This shows 
that even though the states are linearly unstable, they 
influence the flow for a considerable part of its evolution. 
These numbers explain a posteriori why the states could 
be observed in the experiments by Hof et al [lOl . 

2. Upon increasing TZe from 2200 to 2500 the combined 
statistical weight of all detected coherent states decreases 
only weakly. However, there is a clear shift towards states 
with a larger number of streaks. 

3. Due to their prevalence the coherent states sig- 
niflcantly influence the turbulent dynamics at low TZe. 
This opens a route to modelling turbulence by explor- 
ing dynamical interconnections between coherent states. 
To this end we considered the dynamics as a random 
walk between a limited number of coherent states, and 
extracted the transfer probabilities between states from 
the numerical time series. The predictions of the Markov 
dynamics agree very well with the numerically observed 
frequency of occurrence and the life time of the coherent 
states. 

4. The decomposition of the Reynolds stresses into 
contributions arising from irregular motion and contri- 
butions from coherent states with three, four, five and 
six high-speeds streaks allowed us to study the contri- 
bution of different structures to the radial momentum 
transport. Trends in the changes of the radial momen- 
tum transport with TZe could be explained in terms of 
substantial changes of the individual dynamical impor- 
tance (statistical weights) of the states while the proper- 
ties of individual states change only slightly. Both effects 
could be separated based on our statistical analysis. 

We conclude that the methods presented in the present 
paper can be used to quantitatively analyze and describe 
turbulent dynamics close to the transition to turbulence. 
Obviously, they can be extended to projectors providing a 
still more detailed characterization of the flow and used in 
other flows as well. Since the approach does not make use 
of specific features of our numerical setup, it should be 
applicable to the analysis of numerical and experimental 
data alike. 
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